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Up-going and Horizontal Tau Air-Showers, UpTaus and HorTaus, may trace Ultra High Energy Neutrino Tau 
Earth Skimming at the edge of the horizons. Their secondaries (/i* and 7 bundles with pair flashes) might trace 
their nature over UHECR secondaries in horizontal showers. Indeed the atmosphere act as a perfect amplifier 
as well as a filter for showers: down-ward and horizontal /i bundles may still be originated by far Ultra High 
Energy Cosmic Rays skimming the terrestrial atmosphere but their rich gamma component will be exponentially 
suppressed. At large zenith angles after crossing a large slant depth {Xmax > 3 x 10^ g cm~^) the number of fi^ 
and secondary 7's (produced by the pair from fi decay in fiight) is comparable. On the other hand, up-ward 
muon bundles from UpTaus and HorTaus may arise within a young shower with a larger gamma-muon ratio 
(~ 10^), leaving its characteristic imprint. We estimate the UpTaus and HorTaus rate from the Earth and we 
evaluate the consequent event rate of /x*, and 7 bundles. We show that such events even for minimal GZK 
neutrino fluxes could be detected by scintillator arrays placed on mountains at 1 — 5 km and pointing to the 
horizon. The required array areas are within tens-hundreds of square meters. An optimal structure is an array of 
crown-like twin detectors facing the horizons. 

We argue that such detectors will be able to detect both muonic bundles at a minimal average flux of 10~^^ 
cm~^ s~^ sr~^ and electromagnetic particles (7, e*) at 3 x 10~^ cm~^ s~^ sr^^, a few times each year, even for 
the minimal GZK f flux. 



1. Introduction 

The study of ultrahigh energy upward and hor- 
izontal T air showers produced by r neutrino in- 
teractions within the Earth crust has been consid- 
ered in recent years as an alternative way to de- 
tect high energy neutrinos. The i^r UHE sources 
are more common (or Vg) whose tiny masses 
{Arriij > 5 • 10~^ eV) and flavour mixing would 
lead to a i^r appearence. The results are often 
difficult to interpret. The problem of r neutrinos 
crossing the Earth is indeed quite complicated be- 
cause of the complex terrestrial neutrino opacity 
at different energies and angles of arrival. In ad- 
dition, several factors have to be taken into ac- 
count, such as the amount of energy transferred 
in the Vt - t lepton conversion, as well as the r 
energy losses and interaction lengths at different 
energies and materials. This makes the estimate 
of the links between the input neutrino - output 
T air shower very difficult. Such a prediction is 
further complicated by the existence of a long list 
of theoretical models for the incoming neutrino 



fluxes (GZK neutrinos, Z-burst model flux, E^"^ 
flat spectra, AGN neutrinos, topological defects). 
Many authors have investigated this Vr signature 
by their convolution of flux models and the Earth 
new opacity; however the results are varied, often 
in contradiction among themselves, and the ex- 
pected rates may range over a few order of mag- 
nitude (Fargion, Aiello, & Conversano 1999; Far- 
gion 2002; Bertou et al. 2002; Feng et al. 2002; 
Bottai & Giurgola 2003, hereafter BG03; Tseng 
et al. 2003; Bugaev, MontaruH, & Sokalski 2004; 
Fargion et al. 2004; Jones et al. 2004; Yoshida 
et al. 2004). So far, the majority of the current 
studies on this topic is based on Monte-Carlo sim- 
ulations assuming a particular model of the in- 
coming neutrino flux. Most of the authors focus 
on the UpTaus tracks in underground detectors. 

To face such a complex problem, we considered 
the simplest approach: flrst we disentangle the 
incoming neutrino flux from the consequent r air- 
shower physics; therefore, to establish the r pro- 
duction rate we introduce an effective volume and 
mass for Earth-skimming t's, which is indepen- 
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Figure 1. Horizontal Upward Tau Air-Shower 
(HorTauS) originated by UHE neutrino skimming 
the Earth: fan-shaped jets arise because of the 
geo-magnetic bending of charged particles at high 
quota ('^ 23 — 40 km). The shower signature may 
be observable by EUSO just above the horizon. 
Because of the Earth opacity most of the UpTau 
events at angles > 45 — 50° will not be observ- 
able, since they will not be contained within its 
current field of view (FOV). 



dent on any incoming neutrino flux model. This 
volume describes a strip within the Earth crust 
where neutrino/antineutrino-nucleon, i'r{i^T)~N, 
interactions may produce emerging , t+ lep- 
tons which then shower in the atmosphere outside 
the Earth. 

We present a very simple analytical and nu- 
merical derivation (as well as its more sophisti- 
cated extensions) which takes into account, for 
any incoming angle, the main processes related 
to the neutrinos and t leptons propagation and 
the T energy losses within the Earth crust. Our 
numerical results are constrained by upper and 
lower bounds derived in simple approximations 
(see Fargion et al. 2004 for details). The effec- 
tive volumes and masses will be more severely 
reduced at high energy because we are interested 
in the successful development of the r air-shower. 
Therefore we included as a further constraint the 
role of the air dilution at high altitude, where t 
decay and the consequent air-shower may (or may 
not) take place. 

We showed (Fargion et al. 2004) that our re- 
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Figure 2. Column depth as a function of the 
incoming angle having assumed the multi-layers 
structure given by the Earth Preliminary Model 
(see Fargion et al.2004). 



suits give an estimate of the r air-shower event 
rates that exceeds earliest studies but they were 
comparable or even below more recent predictions 
(Yoshida et al. 2004). Secondly, we point out 
here, for the first time, that the consequent /x^, 
e^, 7 signature of HorTaus largely differs from 
that of horizontal UHECR backgrounds. Finally 
we introduce in the calculation of the number of 
events an additional suppression factor related to 
the altitude at which air showers are observed. 
This guarantees the optimal extension and the 
largest flux for the shower to be detected at each 
observational h altitude. 

2. Effective volume and masses for Hor- 
Taus 

To calculate the effective volume we assume 
that the neutrino traversing the Earth is trans- 
formed in a tau lepton at a depth x, after having 
travelled for a distance {D{9) — x). The column 
depth D{9) defined as / p{r)dl, the integral of the 
density p{r) of the Earth along the neutrino path 
at a given angle 9 is shown in Fig. |21 The angle 9 
is included between the neutrino arrival direction 
and the tangent plane to the earth at the observer 
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Figure 3. Comparison between ir(,9) ^nd Ir for 
rock and water. As one can see from the picture, 
iT(^) is shorter than Ir at energies above 10^^ eV, 
thus it corresponds to a smaller effective volume 
where r's are produced while keeping most of the 
primary neutrino energy. The energy label on the 
X axis refers to the newly born tau for ir(/3)- 



location {6 = 0° corresponds to a beam of neu- 
trinos tangential to the earth's surface) and it is 
complementary to the nadir angle at the same 
location. The probability for the neutrino with 
energy E^, to survive until a distance {D{9) — x) 



IS e 



-{D(e)-x)/L, 



, while the probability for the tau 
to exit the Earth is e~^/'^. On the other hand, 
as we will show in the next section, the probabil- 
ity for the outcoming t to emerge from the Earth 
keeping its primary energy is e~'^^^^^'^'> (where 
e-^/LrW) ^ p-^/i. at energy > 3 x lO" cV). 
By the interaction length we mean the char- 
acteristic length for neutrino interaction; as we 
know its value may be associated to the inverse 
of the total cross-section axot = crcc + f JVC, in- 
cluding both charged and neutral current interac- 
tions. It is possible to show that using the acc in 
the e"^^^^-'"^-'/^" factor includes most of the Pr 
regeneration along the neutrino trajectory mak- 
ing simpler the mathematical approach. 

Indeed the use of the total cross-section in the 
opacity factor above must be corrected by the 
multi-scattering events (a neutral current inter- 



action first followed by a charged current one 
later); these additional relay events ("regener- 
ated taus") arc summarized by the less sup- 
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pressing acc factor in the e 
ity term. The only difference between the real 
case and our very accurate approximation is that 
we are neglecting a marginal energy degrada- 
tion (by a factor 0.8) for only those "regenerated 
taus" which experienced a previous neutral cur- 
rent scattering. We also neglected all the charged 
current events deep inside the Earth whose tau 
lepton birth and propagation is source of lower 
energy tau neutrinos which pile up at energy 10^^ 
eV and below ; such a i^r regeneration (absent in 
fij, or i/g cases) may be here neglected because of 
its marginal role in the range of energy (> 10^^ 
eV) we are interested in. Therefore our estimates 
give a lower bound for the detection of t air show- 
ers from UHE Vr skimming the Earth. The rarer 
short range decay of UHE r inside the Earth, 
with a marginal energy loss, may be an additional 
source of "UHE regenerated" Vr even at energy 
10^^-10^^ eV; their additional contribution (here 
neglected) will bo discussed elsewhere. 

The effective volume per unit surface is given 

by 



A ~ 2'kRI 

D(e) 
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where A is any arbitrary surface above the cor- 
responding effective volume. For instance this ex- 
pression has been first estimated for all the Earth. 
In this case A is just half of the terrestrial surface, 
due to the request of selecting only the upward 
direction. 

Under the assumption that the x depth is inde- 
pendent of Li, and It, the above integral becomes: 
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Figure 4. Effective mass for UpTaus and Hor- 
Taus per km square unit area including the sup- 
pression factor due to the finite extension of the 
Earth's atmosphere (an horizontal length of 600 
km). The curves obtained (the red line for an 
Earth outer layer made of water and the pink 
line for the rock) are compared with a simplified 
model of the Earth, considered as an homoge- 
neous sphere of water (green line) and rock (blue 
line). Here we used the interaction length Lr(/3) 
and the volume is expressed as a function of the 
final tau energy. Note that above 10~^ km^ the 
effective mass-volume in the energy range 3 x 
10^^ - 10^^ eV is larger than the volume of the 
atmospheric layer, whose ability to convert down- 
ward neutrino in observable air shower is negligi- 
ble. Only horizontal neutrino interaction in the 
atmosphere may be detected, at a much lower 
rate than the HorTaus ones. 



where the energy of the neutrino E^, has been 
expressed as a function of Et via the introduc- 
tion of the parameter rj = E,y/Erf, the frac- 
tion of energy transferred from the neutrino to 
the lepton. At energies greater than 10^^ eV, 
when all mechanisms of energy loss are neglected, 
rj = Ev/Etj = E^/En — 1.2, meaning that the 
80% of the energy of the incoming neutrino is 
transferred to the newly born r after the v — N 
scattering (Gandhi 1996, 1998). 

When the energy losses are taken into account, 
the final t energy E^^ is a fraction of the one at 



its birth, E^. Their ratio Xi = Etj/Et^ is related 
to 77 by the following expression 
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This ratio defines the fraction of energy of the 
outgoing T compared to the incoming Vt- 

Once the effective volume is found, we intro- 
duce an effective mass defined as 
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where Pout is the density of the outer layer of 
the Earth crust: po„t = 1.02 (water) and 2.65 
(rock). 

Now we can neglect the severe energy degra- 
dation (r/ ^ 1) of the HorTaus, considering only 
those UHE Vr that are converted into r within a 
smaller length ir(/3) i-^T(i3) ^ (c7rir)~"'^ < It) 
(see Fargion et al. 2004 for detail). Note that /3 
is the coefficient due to nuclear photoproduction 
and bremsstrahlung energy losses. 

n(B) D(e) 

Given that in general e e ^.^cc , the 

second exponential in the integral in Eq. [21will be 
also neglected in writing the following equations. 

The expression of the effective volume in this 
most general case becomes 
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where the interaction length ir(/3) (shown in 
Fig. 13 and compared to It) guarantees a high en- 
ergy outcoming t even if outcoming from a thin- 
ner Earth crust (see Fargion et al. 2004 for a 
more detailed discussion of L^(^)). 

The terrestrial chord, D{6), shown in Fig. 12 is 
responsible for the i/j- opacity, and L^, is the inter- 
action length for the incoming neutrino in a water 
equivalent density, where L^^^ = {crcc i^)~^ ■ It 
should be kept in mind that both and D{9), 
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converted into the water equivalent chord, depend 
on the number density n (and the relative matter 
density of the inner shell tracks) . In reality the 
Earth interior has been idealized as an homoge- 
neous sphere of water of column depth D{9) (see 
Fig. [5J and the i^T-interactions are considered for 
the water density. 

We remind that the total neutrino cross section 
CTj/ consists of two main component, the charged 
current and neutral current terms, but the r pro- 
duction depends only on the dominant charged 
current whose role will appear later in the event 
rate number estimate. The interaction lengths 
Lt0, L^ccj depends on the energy, but one should 
be careful on the energy meaning. Here we con- 
sider an incoming neutrino with energy E^^ , a 
prompt r with an energy i?,-. at its birth place, 
and a final outgoing r escaping from the Earth 
with energy , after some energy losses inside 
the crust. The final r shower energy, which is the 
only observable quantity, is nearly corresponding 
to the latter value Etj because of the negligible r 
energy losses in air. However we must be able to 
infer E^-^ and the primary neutrino energy, E^, , to 
perform our calculation. The effective volume re- 
sulting from Eq. ^calculated for a detector with 
a 1 km^ acceptance area is displayed in Fig. 0] 

3. Event Rate of Tau air showers for GZK 
neutrinos with EUSO and Auger 

After having introduced the effective volume 
we can estimate the outcoming event number rate 
for EUSO for any given neutrino flux. The con- 
sequent event rate for incoming neutrino fluxes 
may be easily derived by: 

dA,„ ( f dN, \ 

(5) 

where L^^cc = {crNun)~^, = = 5 • 

10~^* ( lofo'ev ) "^^cm^^ sec^^ sr~^ describes as 
a flat spectrum (a — 2) most of the GZK neutrino 
flux and as a linearly increasing spectrum (a = 
1) the Z-burst model; pr is the density of the 



most external layer (either rock or water). The 
assumption on the flux may be changed at will 
and the event number will scale linearly according 
to the model. 

In Fig. Owe show the expected number of event 
for EUSO where we have included the Earth's 
atmosphere and we have used L^i^p), so that we 
may express the results as a function of the final 
T lepton energy. 

As one can see from Fig. |S1 at energy E — 10^^ 
eV the general expected event rate is given by: 

Ae„ = 5-10-1 W-^s-^sr-i (YmPL] {2tt r]Euso^t) 

V ^i^cc J 

^ V50eycm-2s-isr-i ) ^ eV ) 

where rjEuso is the duty cycle fraction of EUSO, 
rjEuso — 10%, At ~ 3 years. Such number of 
events greatly exceed previous results by at least 
two orders of magnitude (Bottai et al. 2003) but 
are comparable, but below more recent estimates 
(Yoshida et al. 2004). 

We have also calculated the number of event 
rate for the Auger detector (Fargion et al. 2004). 
The number of events is slightly above the unity 
at = 10^9 eV (see Fig. 19 in Fargion et al. 2004) 
for scintillator detection and near 0.1 events/yr 
for photoluminescent detectors. However at en- 
ergy as low a.s E ^ 10^^ eV in GZK model the 
number of events possibly detectable by the scin- 
tillators of Auger increases to 26.3 for a rock 
layer (for scintillator detectors) and 2.6 events 
for photofluorecence ones: the most remarkable 
signature will be a strong azimuthal asymme- 
try (East- West) toward the high Andes mountain 
chain. The Andes shield UHECR (toward West) 
suppressing their horizontal flux. These horizon- 
tal showers originated in the atmosphere at hun- 
dreds of km, are muon-rich as well as poor in 
electron pair and gammas and are characterized 
by a short time of arrival. On the other hand 
the presence of the mountain ranges at 50 — 100 
km from each of the Auger detector will because 
of the mountain geometry by a numerical fac- 
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Figure 5. Number of EUSO Event for HorTaus 
in 3 years record as a function of the outgoing 
lepton tau (-^r(/3) as interaction length), including 
the finite extension of the horizontal atmospheric 
layer. At energy Er = 10^^ eV, the event number 
is Nev = 3.0 {(t>^E^/50 eV cm'^ s"! sr"!) for 
the water and N^v = 6.0 {4i^E^/50 eV cm~^ s^^ 
sr^^) for the rock. The resulting number of events 
has been calculated for an initial GZK neutrino 
flux: cx E-^ . 



tor 2 — 6 (for each given geographical configu- 
ration) the event number respect to the East di- 
rection. These EeVs rare nearby (few or tens of 
km) tau air-shower will be originated by Horizon- 
tal UHE neutrinos interacting inside the Andes. 
These showers will have a large electron pair, and 
gamma component greater than the muonic one, 
and a dilution in the time of arrival (microseconds 
over nanoseconds in horizontal /i^ by UHECR at 
the horizon). Timing, spectra, muon-gamma ra- 
tion and directionality will clearly mark their 
nature. 



4. The differential rate of Tau air showers 

We introduce now a differential expression of 
the number of events which allows to calculate 
the number of events as a function of the angle 
6. We can rewrite the expression of the effective 
volume given in Eq. 0|as a differential volume for 



T T^^ 'K^ -K^ -K^ Vj 



2ir- 



♦ 27 
-> 47 
■+67 



~ 17.4% 
- 17.8% 
- 11.8% 

- 25.8% 
- 10.79% 
- 1.23% 
~ 10% 

- 5.18% 



Unobservable 
1 Electromagnetic 
1 Hadronic 
1 Hadronic, 2 Electromagnetic 
1 Hadronic, 4 Electromagnetic 
I Hadronic, 6 Electromagnetic 

3 Hadronic 
3 Hadronic, 2 Electromagnetic 



Figure 6. Tau decay channels. The probabilities 
for pion secondaries take into account the contri- 
bution of K mesons originated from r decay |S] 



each arrival angle 9: 
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and we obtain the following expression for the 
differential rate of events 
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with 77 = 1.2 and E^„ = lO^^ eV. 



If we now integrate on the solid angle dQ (half 
side) we obtain the above formula multiplied by 
a factor 27r. 

Given the t number of events we can calculate 
the rates of fi, pairs and 7, which originates as 
secondary particles from the t decay (see Fig. EJ. 
The number of muons is related to the total num- 
ber of decaying pions and according to Matthews 
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Figure 7. The differential number of event rate of 
T leptons (HorTaus) for an input GZK neutrino 
flux. As in previous Figures we arc assuming that 
r's are escaping from an Earth outer layer made 
of rock. Note the discontinuity at ~ 1 rad, 
due to the corresponding inner terrestrial higher 
density core (see Fig. ISJ. 
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Figure 8. The differential number of event rate 
of the secondary muons produced by the decay 
in flight of T leptons in the Earth's atmosphere: 
HorTaus. As in previous Figures we are assum- 
ing an input GZK neutrino flux and that r's are 
escaping from an Earth outer layer made of rock. 



(200f ) is given by 



iV„ ~ 3 • f 0^ 
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and we obtain finally 
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We show in Fig. |Slthe average differential rate 
of r's and the secondary fi^, and 7 bundles 
from the decay of r leptons. We find that the 
muon signal at the horizon, related to Earth skim- 
ming tau neutrinos is above 10^^^ - 10^^^ cm^^ 



s~^ sr^-'^. One should notice that the muonic 
background produced by atmospheric neutrinos 
(CR — ^ i^at77i ^ A*^) below the horizon ap- 

proaches at beast to the value — 2 • 10~^^ 

cm~^ sr~^, which is at least one order of mag- 
nitude lower than what we have obtained from 
our calculation for minimal GZK fluxes (see 

Fig. El). 

Moreover this muonic shower would have a sig- 
nificant 7 component, with a high number of pho- 
tons - N^/Nf, - 10^ (Cillis & Sciutto 2001; see 
also Fig. 15 in Cronin 2004) - because of the r de- 
cay channels into both charged and neutral pions 
as one can see from Fig. (H] 

On the other hand, horizontal UHECRs will 
not represent a source of contamination for our 
signal because the horizontal 7's produced in 
the hadronic shower inside the atmosphere would 
be exponentially suppressed at large slant depth 
{Xmax > 3000 g cm^^) and large zenith angles 
{0 > 70°) (Cillis & Sciutto 2001). Only /i± 
can survive when propagating through the atmo- 
sphere at large zenith angles. Such muons would 
also be source of parasite 7 signal - due to the 



8 




Figure 9. The differential number of event rate 
of secondary electron pairs originated from the 
decay of r leptons for an input GZK neutrino flux. 
As in previous Figures we are assuming that r's 
are escaping from an Earth outer layer made of 
rock. 
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Figure 10. The differential number of event rate 
of secondary 7's originated from the decay of t 
leptons for an input GZK neutrino flux. As in 
previous Figures we are assuming that t's are es- 
caping from an Earth outer layer made of rock. 



duced as 



pair produced in the decay in flight - but 
the gamma-to-muon-number ratio would be now 
approximately <^ 1. 

Therefore, this difference would allow to dis- 
tinguish gamma-rich HorTaus from common hori- 
zonal gamma-poor UHECR events to reveal UHE 
earth-skimming i^r^s . 

For a more precise approach to the calculation 
of the rate of events one has also to take into 
account that the number of events varies as a 
function of the height h at which the observer 
is located to detect the muonic, electronic and 
gamma shower. Therefore we introduce an addi- 
tional factor 



dN„ 
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where the parameter hg that we have intro 



hs = Rr{Er) sin0 + sing (12) 

Pr 

defines the optimal height where the shower can 
reach its maximal extension at the correspond- 
ing energy Er- This is the sum of the height 
reached by the r in the atmosphere before its 
decay (_R^sin0), where we have neglected t en- 
ergy losses in the atmosphere, and the altitude 
reached by the secondary particles of the shower 
which is related to the parameter Xmax- Note 
that Rr = A.9{Er/lO'^'^ eV) km and X^ax/pr = 
5.7 -I- 0A6\n{Er/10^'^eV) km. Here we have con- 
sidered the air density pr = 1.25 x 10"'^ g cm~^ 
constant and equal to the value at the sea level. 
At low energies (10^^ - 10^^ eV) the second term 
is dominant and pr can be considered as constant 
because the r lepton travels for less than 1 km be- 
fore it decays. At higher energy (10^^ - 10^^ eV) 
the first term is dominant, and we can neglect the 
way the exact value of pr changes with the alti- 
tude. We show our results in Fig. II II and Fig. IT^ 
where the secondary fi^ fluxes are described for 
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Figure 11. The average differential event rate of 
secondary muon leptons from HorTaus having as- 
sumed Et = 10^^ eV and an input GZK-hke neu- 
trino flux at different altitudes. Again we are as- 
suming that t's are escaping from an Earth outer 
layer made of rock. Note the discontinuity of the 
angular spectrum due to the inner core of the 
Earth at ~ 1 rad and the asymptotic behaviour 
for ft, 1 km. 
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Figure 12. The differential event rate of muon 
leptons having assumed Er = 10^^ eV and an 
input GZK-like neutrino flux. Again we are as- 
suming that t's are escaping from an Earth outer 
layer made of rock. Note that because the Earth 
is nearly opaque to EeV , only Earth-skimming 
neutrinos nearly horizontally are visible. The dis- 
continuity in the spectrum at ~ 1 rad occurs 
also at energies as large as 10^^ eV, but at lowest 
flux, and we do not see it in this picture. 



two given Er energies. 

5. Conclusions 

Horizontal showers from normal hadrons (or 
gammas) are strongly depleted of their electro- 
magnetic component because of the large slant 
depth {Xjnax > 3 X 10"* g cm~^), while horizontal 
tau air showers are not. Indeed "young" Hor- 
Taus either of hadronic (67%) or electromagnetic 
(33%) nature (see Fig. at their peak shower ac- 
tivity are expected to have a large N^/N^ ratio, 
greater than 10^ (but with a characteristic en- 
ergy ratio E^/Efj_ < 10^). Old horizontal showers 
would have N^/N^ ~ 1. This difference would al- 
low to distinguish and disentangle HorTaus from 
horizontal UHECR events, even in absence of 
good angular resolution, opening a new perspec- 
tive in the UHE neutrino Astronomy. The sec- 



ondary fluxes of muons and gamma bundles made 
by incoming GZK neutrino fluxes and their Hor- 
Tau showers, is well above the noise (by one-two 
order of magnitude) made by up-going muons, 
Earth-Skimming trace of atmospheric neutrinos. 
The neutrino signals at energies much above EeV 
may be even better probing the expected harder 
neutrino Z-Burst model spectra [H], PT] . 

The peak fluences we find in the /i and 7 com- 
ponent at the horizon (±5°) will give a signal 
well above the background produced by atmo- 
spheric J^'s. We did not discussed the albedo 
muons whose fluxes measured by Nemo-Decor ex- 
periments, (j)^ <^ 10~^ cm~^ s~^ sr~^, are made 
mostly by single tracks ^HI- Because pair (or 
triple) bundle muons are much rarer {(piipaiT < 
10^^(/)^^.,,g,^ ~ 10^^'^ cm~^ s~^ sr~^), the search 
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and detection of muon bundles by GZK HorTaus 
at a minimal rate of 10^^^ cni^^ s^^ sr~^ (over 
an area of 10^ m^) will lead, in a year, to about 
30 muons possibly clustered in five - ten multiple 
bundles. These events will be reinforced by hun- 
dreds or thousands of associated coUinear gamma 
flashes. A detector with an area of few tens or 
hundreds of square meters pointing to the hori- 
zon from the top of a mountain would be able 
to reveal the GZK — t young showers (fori, 
Sergi & Fargion 2004). The characteristics of a 
prototype twin crown-like array detector to be 
placed on mountains, balloons, or satellites ^Hl 
will be discussed in detail elsewhere. The simul- 
taneous sharp 7 bundle at (p^ ^ 10~^ 10~^^ 
cm~^ s"^ sr^^ and the "burst" of electron pair 
at 4)e+e- ^ 10^^ cm^^ s^^ sr^^ would give evi- 
dence of unequivocal r signature. It should be re- 
mind that the neutrino interaction enhancement 
by TeV new Physics (as shown in FigT and Ap- 
pendix C) would produce also an increase of 
hundreds or thousands time in HorTaus beyond 
a mountain Chain (like Auger) than standard 
weak interactions would do. Therefore Auger 
must soon detect either the Ande Shadows for old 
UHECR (at zenith angle larger than 85° - 88°) 
by their absence as well as young [7] , P HorTaus 
productions by the mountain themselves due to 
New Physics at TeV. In a few years Auger might 
be even able to observe also GZK neutrinos in- 
duced HorTaus at EeV energies (by sure if some 
technical improvement will be made). EUSO ex- 
periment will be able to see at least half a dozen 
of events of HorTaus mainly enhanced along the 
Continental Shelves or Mountain edges. To con- 
clude we want to remind that inclined-vertical 
PeVs r air showers (UpTaus) would nearly always 
be source of 7 "burst" surviving the atmosphere 
opacity. These sharp UpTaus (with their com- 
panion HorTaus above and near EeVs) might be 
observed by satellites as brief Terrestrial Gamma 
Flashes (TGF). Indeed we identified a possible 
trace of such events in BATSE record (taken dur- 
ing the last decade) of 78 upgoing TGF maybe 
associated with galactic and extragalactic UHE 
neutrino sources [H|- 
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